Line of Dirac monopoles embedded in a Bose-Einstein condensate 
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The gauge field of a uniform line of magnetic monopoles is created using a single Laguerre-Gauss 
laser mode and a gradient in the physical magnetic field. We study the effect of these monopoles 
on a Bose condensed atomic gas, whose vortex structure transforms when more than six monopoles 
are trapped within the cloud. Finally, we study this transition with the collective modes. 
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A point source of magnetic flux, a magnetic monopole, 
has long been an important missing component of grand 
unified and superstring theories. Despite an exhaustive 
search for such monopoles [11 12] , no convincing evidence 
for the existence of a stable magnetic monopole has been 
forthcoming. Analogies to certain aspects of monopoles 
have been uncovered in physically tractable systems in- 
cluding spin ice [3 , topological insulators [4 , the anoma- 
lous quantum Hall effect [5 , and superfluid ^He [6]. Nev- 
ertheless, the experimental realization of a true Dirac 
monopole remains an important target. 

Over the past three years the cold atom gas has 
emerged as a promising new system in which to explore 
the motion of particles in magnetic fields |7]. Starting 
with a seminal experiment [8, 9 that imposed an effec- 
tive uniform magnetic field on a Bose-Einstein conden- 
sate (BEG) of neutral atoms, the experimental protocol 
has been extended to realize spin-orbit coupling [10 , and 
effective electric fields [11 . Following these successes, the 
spotlight has recently turned to generating the high mag- 
netic fields required to study strongly correlated states 
such as the fractional quantum Hall effect p!2HT8] . How- 
ever, cold atom gases also present a unique opportunity 
to study magnetic field configurations not realizable in 
the solid state, in particular a magnetic monopole [19]- 
[22] . The experiments suggested to date demand either a 
laser setup of a superposition of Laguerre-Gauss modes 
along with Hermite-Gauss modes [T9H2T] . or three mag- 
netic field modes [22], which have not yet been real- 
ized in practice. Moreover, the configuration of a single 
monopole of fixed charge means that there is no tuning 
parameter that would allow the direct characterization 
of the state [2T]. 

Here we propose a significantly simpler experimental 
protocol to embed a uniform line of monopoles into a 
BEG. Motivated by previous experiments [8l[T0l[TTl that 
successfully realized effective gauge fields, our scheme 
adopts the same physical magnetic field gradient, and 
replaces the co-propagating laser plane waves with the 
first optical Laguerre-Gauss mode [23H25] . The system 
undergoes a transition with increasing monopole density, 
and we characterize the two phases by their collective 



modes spectra. 

We first outline the experimental scheme to create the 
gauge field of a line of monopoles. Starting from the ge- 
ometry in Fig. [ija), we trap bosonic atoms that have 
three coupled internal levels (for example the F = 1 
manifold of ^^Rb) labeled by mp. We apply a small 
physical magnetic field {B + B'z)z where the magnetic 
field gradient B' could be generated by an anti-Helmholtz 
coil pair. The field introduces a linear Zeeman split- 
ting so the atomic levels have energies gfiBm^B -\- mp^, 
where 6 = gji^B' z is the detuning from the Raman res- 
onance as shown in Fig. [TJc) . Go-propagating Laguerre- 
Gauss laser beams induce the complex Rabi frequencies 
= K.o{V2p/VIwYexp{£/2 - p'^/w'^ ± Ucj) + \kz) [7], 
with £ the beam winding number, w the beam width, 
p = \/x'^ -\- y'^ the distance from the z-axis, (j) the az- 
imuthal angle, hcq the Rabi frequency at p = w^/i/2^ and 
k the plane wave vector. The single-photon detuning A 
from the excited states manifold should be large versus 
the Rabi frequencies A so that the excited states 

can be adiabatically eliminated and the Hamiltonian acts 
only on the ground state manifold. 

The Hamiltonian in the basis set {\mF = — 1), \mF = 
0), \mF = +1)} is ^ = (-n2vV2m+V'-/i)l+/7, with m 
the atomic mass, the spin-independent trapping potential 
is V = muj'^{p'^ -\-z'^)/2 with uj the trap curvature, and in 
the rotating wave approximation the coupling matrix is 

^ I -25 \ 

[/ = - 6^.0/^+ , (1) 

^ \ 25 ) 

with e the quadratic Zeeman shift in units of kq. The 
three bands of the Hamiltonian are split hy ^ ^ 5. 
If the thermal energy and chemical potential p [26 are 
much lower than hn^ [8 then only the ground state is 
occupied 

lx)=e^^^^cos2 \\ - l)-^|0)+e-2«^sin2 , (2) 

with mixing angle tan^ = \k+\/\/25, and energy 
-h^5^ + kI/2. 
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FIG. 1: (Color online) (a) Geometry of the proposed scheme. 
Co-propagating laser beams with winding numbers =b^ and 
the physical magnetic field B are aligned along z. (b) The 
gauge field realized, with the line of monopoles in orange, 
the effective magnetic field as green arrows, and the atomic 
cloud as the blue sphere, (c) The level diagram for the three 
electronic ground states displaced by the external physical 
magnetic field B. The irradiating laser beams and —I 
induce resonant coupling. In practice A ^ ^. (d) The relative 
axial to radial magnetic field Bp. The dashed green line 
denotes the path of constant radius in the trap, and the blue 
sheet is B^/Bp = 0. 



Following the prescription of Ref. [7] , we can calculate 
the gauge field as A = 2\h{x\^x) = M cos 0(f)/ p. The 
corresponding effective magnetic field is B = V x A so 



MS' 



( 



a+/3i 



(3) 



where a = \/2 + e + 3eVlO and /3 = l/\/2 + e/2^ 
17eVl28. The scalar potential defined hyW = h?[{Ve)^-\ 
sin^ 0(y(p)^]/8m is 




(4) 



After expressing the magnetic field gradient through S = 
S'z^ we expand these expressions in the /^o ^ S'w limit, 
and about the radius Pmax = of maximum laser 

beam intensity to yield the effective magnetic field 




(5) 



and in the same limit the scalar potential is W = 
h^{i/ p)'^ /2m. The effective magnetic field formed has the 
configuration shown in Fig. [ijb) : it is directed radially 
and decays as so it represents the magnetic field em- 
anating from a line of Dirac monopoles along z, with line 
density 2\/2le^5' / kq. However, from Eqn. (|3|, = {) 
only exactly at the radius Pmax = i^\/-^/2, but remains 
small away from it. To minimize B^ over all radii we set 
e = —2.8, which yields the variation of B^/Bp shown in 
Fig.jlTd). Over the range of p and z within a trap (in 



Fig.ji]; 



d) set so p^ < w'^) we find that \Bz/Bp 



< 0.1 



and so the effective field closely approximates that of a 
line of monopoles. We also note that the anti-Helmholtz 
coil pair necessarily generates a radial physical magnetic 
field gradient as well as the B' along the z-axis. This 
generates a Bz component that could also be removed by 
the quadratic Zeeman shift. In both cases as A is con- 
tinuous away from p = flux is conserved and the only 
monopoles lie along the z-axis. The Bz component of the 
effective magnetic field is then a small field component 
corresponding to a changing line density of monopoles. 
This will be further tested in Fig. [2] by comparing the 
BEG ground state for both gauge fields. 

To probe the experimental prospects of the proposed 
configuration we first check the maximum achievable 
monopole density. For ^ = 1, using the realistic [8] 
values S' = A.bE^/pm and /^o = 16£^r, with the 
recoil energy, a cloud of diameter d = 80pm contains 
N = 2^/2eS^d/ hzo ^ 80 monopoles. We adopt the natu- 
ral definition of monopole charge density so that the flux 
quantum emanating from a single charge winds the 
phase of the BEG by a factor of 27r. 

Having formed the effective gauge field of a line of 
monopoles, we now develop the formalism to study its in- 
fluence on a BEG. To seek and study the ground state we 
discretize [211 El] the Hamiltonian onto a lattice of spac- 
ing a. The derivative at a lattice site R = {x, z} along 
a single link r is transformed to [Drtl^(Il)]^[Drtl^(Il)] 

(V^R+r - t^R,rV^R)UV^R+r " ^R,rV^R)^^ C^R^r, whcrC 

the hopping matrix is Un,r = exp(iA(R) • r/fi), and t/^r 
is the condensate order parameter at R. With this in 
place the mean-field energy is 



R,r R 



with external potential Fr = mcj^(x^+?/^ + z^)/2, and a 
chemical potential p to fix the number of particles. We 
minimize the free energy variationally with respect to tp^ 
to yield the discrete Gross-Pitaevskii (GP) equation 



2ma^ 

r 

+ (Vr + W^r + c|VrHVr- (6) 

To calculate the ground state we solve the GP equations 
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FIG. 3: (Color online) Flux $ exiting the cloud as a function 
of the number of embedded monopoles N. The red line (+) 
shows the flux exiting radially, the blue line (x) the flux ex- 
iting along the Dirac string. The grey shading denotes the 
region where the escaping flux is necessarily less than the 
monopolar charge contained. The vertical black line shows 
the theoretically predicted crossover when flux first exits away 
from the string. 



FIG. 2: (Color online) The particle density isosurfaces of 
BECs containing (a) AT = 4, (b) iV = 12, (c) iV = 55 
monopoles for the gauge field of the proposed setup. The 
isosurfaces for the artificial situation of a perfect line of (d) 
= 4, (e) = 12, and (f) AT = 55 monopoles. The bright- 
ness denotes the depth of the site on the 2/- axis. The right- 
hand plots of (a,b) view the condensate along the z-axis, and 
use the hue to display the condensate phase. 

iteratively using the over-relaxation scheme on a grid of 
200 X 200 X 200 sites with periodic boundary conditions. 

With a tool to determine the ground state wave func- 
tion in place, to calculate the collective mode spectrum 
we adopt the linear-response theory [27 that was suc- 
cessfully used to study the excitation frequencies of a 
uniform BEG in a vortex state [28]. We first substitute 
the general form for the solution of the Gross-Pitaevskii 
equation 

*R,t = e'^*/^ [^R + WA.Re-'"^* + «I,Re'"^*] , (7) 

into the Hamiltonian. Here is the ground state wave 
function and A is the index of the eigenvalue Qx with 
eigenfunctions iiA,R and vx^n- This recovers the time- 
independent Gross-Pitaevskii equation for t/jn along with 
the linear-response Bogoliubov equations 

f c,^^,2 -^l^nf"''''')=^^Af"^'^V (8) 

VflV^Rr -^n J \vx,ii J \vx,n J 

with = -h^ Er(^R-r,r ' 2 + Un,r)/2ma^ + + 

— M + cIV^Rp. Using the condensate order parameter 
for the ground state, we use the Arnoldi method to solve 
the Bogoliubov equations to find the lowest eigenstates 
and eigenmodes that are the collective modes. 

With the physical setup and numerical procedure in 
place, we now study the ground state configuration of 
our trapped BEG embedded with a line of monopoles. 
We first examine the phase behavior of the ground state 
with increasing monopole density. Second, to fully char- 
acterize the state, we calculate their collective mode spec- 
trum. We study 100,000 atoms in a spherical trap with 



= 0.1 and interaction strength c = 0.05 in atomic units 
h = m = 1. The Laguerre- Gauss beams have winding 
number i = 1 and beam width w = 40a, and vary the 
detuning gradient S' / f<io to control the monopolar charge 
density, which within a characteristic length w of the 
center is 4:^/2eS^w/ hZQ. We use the realistic gauge field 
A = 2hcos6^/p. 

With no embedded monopoles the BEG is spherically 
symmetric in the trap. In Fig. |2|a) we examine the den- 
sity isosurfaces for a cloud that encloses N = 4 monopo- 
lar charges. The density falls along the z-axis and the 
phase of the condensate order parameter winds by An 
around the z-axis showing that the flux escapes symmet- 
rically through the top and bottom of the cloud forming 
vortices. On increasing the number of trapped monopoles 
to N = 12, in Fig. [2|b) now three flux quanta pass 
through the top of the cloud, and a further three through 
induced vortices directed away from the axis in the up- 
per hemisphere, with the situation repeated in the lower 
hemisphere. The stepwise change where flux starts to 
escape away from the z-axis occurs dit N ^ 6. The 
cylindrical symmetry of the cloud is broken by the ran- 
dom site updates on the computational lattice. When 
a large number of monopoles N = 55 are trapped by 
the cloud in Fig. [2|c) we observe that the vast major- 
ity of the flux escapes through the side of the cloud 
broadly symmetrically. Finally, we compare the predic- 
tions for the gauge field realized by the proposed laser 
setup, A = 2hcos6(p/ p (Fig. |2|a,b,c)), with an idealized 
A = Aeh54>/nQp of the artificial situation of the magnetic 
field of a perfect line of monopoles (Fig. [2|d,e,f)). The 
distribution of the vortices formed by the exiting flux are 
broadly similar. Therefore the gauge field realized by 
the laser setup accurately models that of an ideal line of 
monopoles. We also checked the stability of the results 
while changing the random number seed of the lattice 
updates, boundary conditions between periodic and hard 
wall, the orientation of the line of monopoles relative to 
the lattice, and lattice size. 

Having studied clouds containing different monopolar 
charges, we now summarize their behavior in Fig. |3] In 
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FIG. 4: (Color online) The collective mode frequencies renor- 
malized by the trapping potential ^x/uj as a function of 
number of embedded monopoles N for different interaction 
strengths and trapping potentials. The dashed arrows high- 
light the expected modes in the = and N ^ oo limits. 



particular we focus on the route of the escaping flux, 
which is along the 2:-axis for N < 5 and penetrates 
through the side walls of the cloud for > 8. When 
A/" > 8, the flux still exits along the z-axis, but remains 
fewer than a total of ~ 6 flux quanta. To understand 
why there is a transition in behavior at 5 < A/" < 8 we 
can compare the energy cost of flux exiting the cloud 
through the two different routes. The energy of a vor- 
tex per unit length is [36 where 7 is the number 
of circulation quanta of the vortex and C a constant de- 
fined in [36 . For a system of radius length L, and 
monopole line density 7^, the energy of all the flux es- 
caping along the z-axis is Crf / 12^ whereas if the flux 
exits radially the total energy is CrjLR. Therefore the 
flux exits radially when 12R > rjL^. Here the cost of 
multiple flux quanta exiting along the same vortex line 
out of the top or bottom of the cloud is outweighed by 
the energy cost of flux exiting through longer vortex lines 
out of the cloud sides. For a spherically symmetric trap 
where L = 2R and N = rjL, the flux exits radially when 
A^ > 6. This can be readily seen in Fig. [3] where flux exit- 
ing radially rises rapidly when more than six flux quanta 
are enclosed by the cloud. Here, as in previous studies 
of rotation-induced vortices in BECs, the vortices will be 
too small to image directly. However, if the cloud is bal- 
listically expanded the vortices are magnified and could 
be imaged [29H32 . 

On top of direct imaging, a second experimental tool 
to characterize the BEC with embedded monopoles is 
the collective mode spectrum. This was previously used 
to analyze both the uniform and rotation-induced vor- 
tices in BECs [27H291 \33\ [34]. Collective modes can be 
excited by a red-detuned optical dipole potential that 
draws atoms into the center of the condensate [29 . 

To study the collective mode spectrum of the cloud 
we take the ground state wavefunction and diagonalize 
the linear response Bogoliubov equations. In Fig. |4] we 



show the lowest mode frequencies for a series of clouds 
with differing numbers A^ of embedded monopoles. We 
first compare the observed spectra at A^ = to that 
of a vortex-free BEC in a spherical trap whose collec- 
tive modes in the weakly interacting limit are simply the 
single-particle modes in the parabolic trapping potential, 
so Qx = Xu with A G Z+. These frequencies are in good 
agreement with the calculated modes, and double in en- 
ergy on doubling uj from (b) to (c), but barely change on 
reducing the interaction strength c from (a) to (b). On 
increasing N < 6 the flux quanta initially escape along 
the z-axis, breaking the spherical symmetry, and splitting 
the degenerate modes. With a large number of embedded 
monopoles A^ ^ 00 the surface of the cloud resembles a 
two-dimensional BEC penetrated by a vortex lattice (see 
for example Fig. [2|c)). The modes of such a sheet are 
Q = h[nc{nl^nl^nl)/mY^'^ [35^, with {n^,ny,n^} G Z, 
n the density, and L the sheet size, here the circumfer- 
ence of the cloud. These modes were verified by their 
frequencies halving on reducing the interaction strength 
from (a) to (b). This provides an intuitive picture of the 
low energy collective modes in the N ^ 00 limit where 
the surface modes dominate. With these signatures the 
collective modes could help diagnose the behavior of the 
BEC embedded with monopoles. 

In this paper we have constructed a gauge field of a line 
of magnetic monopoles embedded in a BEC. The BEC 
underwent a transition from channeling the flux along the 
Dirac string to the flux escaping out of the cloud walls. 
Finally, we characterized this transition through the col- 
lective modes spectrum. However, there are opportuni- 
ties to explore monopolar gauge fields in two other sys- 
tems. The first is to study the single atom bound states in 
a non-interacting gas. These feel the same effective mag- 
netic field as the condensate atoms, and could be charac- 
terized by exciting transitions between different states. 
The second possibility is to generalize our single spin 
species BEC up to a two-component BEC by exploiting 
a multi-level system with doubly degenerate dark states 
and study the emergence of topological states [2T] . 
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